Initial-state randomness as a universal source of decoherence 
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We study time evolution of entanglement between two qubits, which are part of a larger system, 
after starting from a random initial product state. We show that, due to randomness in the initial 
product state, entanglement is present only between directly coupled qubits and only for short times. 
Time dependence of the entanglement appears essentially independent of the specific hamiltonian 
used for time evolution and is well reproduced by a parameter-free two-body random matrix model. 
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I. INTRODUCTION 

In the last few years quantum entanglement is one of 
the most active research fields in quantum physics, for 
a review see Refs. On microscopic level of coherent 
quantum systems entanglement can be used in various 
quantum protocols to perform non-classical operations. 
On the other hand, in the world of macroscopic objects 
with many degrees of freedom at high temperatures there 
is apparently no observable manifestation of entangle- 
ment. There have been many attempts to explain this 
classicality of macroscopic systems. Most notably, de- 
coherence due to external degrees of freedom is usually 
credited as being responsible for the disappearance of en- 
tanglement from macroscopic superpositions. In a nut- 
shell, the argument goes as follows: even if the system is 
in an entangled state at the beginning, e.g., in a coherent 
superposition of two macroscopic states, time evolution 
will in general transform this coherent superposition into 
an incoherent (i.e., classical) mixture. The reason for 
such decoherence is an always present residual coupling 
of our central system to many uncontrollable external 
degrees of freedom - the environment. For a review of 
decoherence see Ref. However, one must be aware 
that the evolution of the central system plus environ- 
ment is still unitary and therefore, even though the final 
state of the central system and environment will presum- 
ably be very complex, in principle, it will be a pure state 
possessing some bipartite entanglement. 

The resolution of this apparent paradox is similar to 
the one with the second law of thermodynamics 
Increasing of the thermodynamic entropy with time is 
seemingly in contradiction with the reversibility of the 
underlying equations of motion. For explanation one 
can use two observations: (i) practicality - performing 
time reversal by, e.g., reversing velocities of all particles 
might be close to impossible from a practical point of 
view; (ii) probability - initial conditions are prevalently 
of such form that in almost all cases the entropy will 
subsequently increase. Similar arguments can be used 
to explain the apparent lack of entanglement in quan- 
tum systems with many degrees of freedom. First, even 
though joint pure state representing central system and 
environment is bipartite entangled, the detection of en- 



tanglement might be close to impossible because it would 
require very complex measurements involving very many 
degrees of freedom. Indeed, using entanglement wit- 
nesses it has been shown that the detection of entangle- 
ment in a sufficiently complex state gets exponentially 
hard with increasing number of particles [4j. For all 
practical purposes the detection of entanglement in such 
states is impossible. Second argument, that is the role 
played by initial conditions in the course of loosing en- 
tanglement by time evolution is the subject of present 
work. 

We are going to study how the entanglement between 
two qubits changes during hamiltonian time evolution. 
Hamiltonian evolution will act on a system of totally n 
qubits, two of which will be chosen as our central system 
of interest while the remaining n — 2 will act as the "en- 
vironment" . The idea is to study how a general hamilto- 
nian evolution changes entanglement of a smaller subsys- 
tem, whose degrees of freedom we presumably are able 
to control and therefore also measure its entanglement. 
Time evolution with a general hamiltonian, say quan- 
tum chaotic one, will in general produce states whose 
statistical properties are well described by those of ran- 
dom quantum states. For random quantum state on n 
qubits one knows Q that tracing out n — 2 qubits will, 
for large n, with high probability result in a separable 
two qubit reduced density matrix. Therefore, sufficiently 
"complex" time evolution will eventually wipe out en- 
tanglement between two qubits. How are things then, 
for instance, in integrable systems, which in general do 
not generate completely random states? One point we 
have not touched so far is the role played by initial con- 
ditions. For integrable systems there can exist simple 
initial states for which entanglement will persist also for 
long times, nevertheless, as we will see, the majority of 
initial conditions is such that entanglement between two 
qubits will rapidly decay with time irrespective of the 
hamiltonian. This universality will be a consequence of 
the generic form of initial states - their randomness. 

The initial pure state will be chosen to be a product 
state on the central system (two qubits) and either ran- 
dom or random product state on the rest. Therefore, 
initially there will be no entanglement between the two 
chosen qubits. We arc then going to study how much en- 
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tanglement can be produced by various hamiltonian evo- 
lutions and how long will it take until it disappears. En- 
tanglement will change with time due to two competing 
effects. One is entanglement production due to time evo- 
lution with non-separable hamiltonian, while the other 
is entanglement loss due to the spreading of initial state 
randomness throughout the system and the approach of 
system's state to a random state. The net result will be 
the increase of entanglement at short times and a com- 
plete lack thereof after some critical time. In addition, 
time dependence of entanglement will turn out to be uni- 
versal, that is independent of the specifics of the hamil- 
tonian used in time evolution. 

There have been many studies of entanglement evolu- 
tions, let us here mention only those that are closer to the 
present work and deal with two qubit entanglement. Evo- 
lution of concurrence for initial product states has been 
studied for XY model in magnetic field in Ref. Q , see also 
Ref. 0. For sufficiently strong coupling between qubits 
concurrence initially increased with time after which it 
rapidly decreased to zero, similarly as in the present 
work. Evolution of entanglement for initially entangled 
states of two qubits (Bell states) which are weakly cou- 
pled to a generic environment has been studied in Q. 
In such cases entanglement monotonically decreases with 
time from its maximal value at time zero. Evolution of 
entanglement for initial Bell state in an integrable XY 
model has been studied in 0. Our hamiltonian will be 
homogeneous in space and therefore the two qubits of the 
central system will be coupled. This must be contrasted 
to studies of the so-called environment induced entangle- 
ment generation [icj |. where two qubits are uncoupled. 
Disappearance of entanglement (on average) after finite 
critical time found in the present work is reminiscent of 
the so-called sudden death of entanglement, where ini- 
tially entangled state of two uncoupled qubits becomes 
separable after a finite time of open system dynamics [ill ] . 
A system consisting of two spins coupled to electrons has 
been studied for initial separable states in Ref. [T3] , while 
time evolution of entanglement for initial thermal states 
in a XY model has been considered in Refs. fl3l. UM . For 
concurrence in a kicked Ising model see [lj| . In [l6j it has 
been found that starting from an initial non-disordered 
product state in a spin glass model entanglement can 
persist for long times. Von Neumann entropy of a block 
of spins in Ising model and for product initial state has 
been studied in fl7ll. 



II. QUANTIFYING ENTANGLEMENT 

We arc going to study entanglement between two 
qubits which are in turn part of a larger n qubit sys- 
tem. For two qubits positivity of partially transposed 
density matrix with respect to one qubit, p TA , is a nec- 
essary and sufficient condition for its separability [l8l ]. 
Negative eigenvalues (for two qubit states there can be 
at most one) therefore signal the presence of entangle- 



ment. A quantity measuring this is negativity (l9l | N(p) 
which is equal to the sum of absolute values of negative 
eigenvalues of p TA and can be defined as 



N(p) 



(1) 



with the trace norm \\A\U = tr V At A. For two qubit 
density matrices it is simply N(p) = |A^ n | if the minimal 
eigenvalue A^ n is negative and otherwise. 

Entanglement of formation (20| , which quantifies quan- 
tum resources needed to create a given state, has nicer 
mathematical properties as negativity (or logarithmic 
negativity). For two qubit systems entanglement of for- 
mation E-p(p) can be calculated in terms of a simpler 
quantity called concurrence C(p), defined as 



C(p) = max{0, Ai — A 2 — A3 — A4}, 



(2) 



where A; are square roots of decreasingly ordered eigen- 
values of (pa v (g> a v p*a y ® a v ), calculated in standard 
computation basis. Ep(p) is then given by 21 1 



Ef(p) = H 



1 + V /1-C 2 ( P ) 



(3) 



with H(x) = — cclog 2 x — (1 — x) log 2 (1 — x) being a bi- 
nary entropy. For pure states the entanglement of forma- 
tion is given by the von Neumann entropy of the reduced 
density matrix, while it is defined by a convex roof ex- 
tension (minimization over all convex realizations of p) 
for mixed states. A state is separable iff its concurrence 
or iff its negativity is zero. 

The third and last quantity used in measuring entan- 
glement will be the so-called fully entangled fraction [2(| 
defined as 



f(p) = max(4>\p\ip), 



(4) 



where maximization runs over all maximally entangled 
states obtained by local unitary transformations from 
maximally entangled state, i.e., = U\ ® E/2GOO) + 
\11))/V2. One can distill maximally entangled singlets 
from an ensemble of p using BBPSSW [HJ protocol iff 
/ > 1/2. Fully entangled fraction / can be used as a 
lower bound for the entanglement of formation [20| . 

E F > h(f), (5) 

where h(f) is expressed in terms of binary entropy H(x), 



Kf) = 
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In the above inequality an equal sign holds if p is pure 
state. Note that / < 1/2 does not necessarily mean 
that the entanglement of formation is zero. Fully en- 
tangled fraction also determines maximal teleportation 
fidelity [23|. Fully entangled fraction f(p) is equal to the 
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largest eigenvalue of the real part of the density matrix p 
written in the Bell basis, in which all maximally entan- 
gled states have real expansion coefficients. If / > 1/2 
this is in turn equal to [24| 

f(p) = \(i + e( P )), e(p) = ||r|| 1 = trVrtr > (6) 

where T is a real 3x3 dimensional correlation matrix 
given by T,j = tr (p a 1 ® a J ) , with a 1 being Pauli matri- 
ces, i.e., i,j G {x,y,z}. Because of its simple analytical 
form we are going to study 6(p) © rather than f(p) 
Q. They essentially give the same information in the 
interesting regime of / > 1/2. If O(p) > 1 then the state 
p can be used in entanglement distillation. Norm of the 
correlation matrix T can be used as a simple entangle- 
ment criterion also for many-qubit systems (25[. 

In the following we are therefore going to study neg- 
ativity N(p) (sometimes just minimal eigenvalue of the 
partially transposed matrix A^f n ), concurrence C(p) and 
Q(p) which is connected to the fully entangled fraction. 
As we will see, qualitatively all behave in the same way. 
From the analytical viewpoint though O(p) is the sim- 
plest quantity and is therefore the best candidate for an 
analytical treatment. 

III. SYSTEMS STUDIED 

We are going to study entanglement evolution for var- 
ious one dimensional spin hamiltonians consisting of n 
spin-1/2 particles. To verify that the results do not de- 
pend on the underlying dynamics we will use both chaotic 
and intcgrablc systems. 

Heisenberg spin-1/2 model is an intcgrablc model with 
hamiltonian 

H = Y,crfaf +1 + ^ay +1 +a?a? +1 . (7) 

i 

We have checked that the results are similar for 
anisotropic Heisenberg model as well as for isotropic 
Heisenberg model in a tilted magnetic field. 

We can break integrability of the Heisenberg model by 
applying magnetic field, for instance, a staggered field in 
z-direction, 

H = E <<f i + + + E (8) 

i i 

where the strength of the magnetic field is h,2% = and 
/i2j+i = ~\ on odd sites. As we can see in Fig. [1] spac- 
ing of neighboring energy levels agrees with the Wigner- 
Dyson distribution, p(s) = S7r/2exp (— s 2 tt/A), which ap- 
proximates distribution of spacings for gaussian orthogo- 
nal random matrix ensemble and is typical for quantum 
chaotic systems [26l |. 

Last model will be Ising chain in tilted magnetic field, 

// E«: ' < ' ( 9 ) 
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FIG. 1: Level spacing distribution for the Heisenberg model in 
a staggered field (O for n = 16. Dashed line is Wigner-Dyson 
distribution holding for quantum chaotic systems. 

Ising model in tilted magnetic field also displays typical 
signatures of quantum chaos [27| . similarly as its time- 
dependent kicked version [28[ . We have checked that the 
results are essentially the same also for intcgrablc trans- 
verse Ising model. 

As we will see, all three models will display similar 
evolution of entanglement being in turn equal to the one 
for a two-body random matrix model. Therefore, our 
main focus will actually be on a two-body random matrix 
model, for which we have only nearest neighbor coupling 
terms, 

H = Y J Ki+^ ( 10 ) 

i 

with hij+i acting nontrivially only on two qubits, for 
which it is a 4 x 4 random hcrmitian matrix, same for all 
coupled pairs and normalized as tr(ft,f i+1 ) = 1. A ran- 
dom hermitian matrix is a matrix whose matrix elements 
are independent random complex gaussian numbers [2(| . 
We always averaged over an ensemble of random matrices 

For all hamiltonians the geometry is that of an one- 
dimensional chain with open boundary conditions. The 
state at time t is obtained as \4>{t)) = exp (—iHt)\ip(0)) 
from which we get the reduced density matrix for the two 
qubits between which we study entanglement, 

P (t) = tr n _ 2 |V(t)>W*)l. (11) 

where a subscript n — 2 means tracing over n — 2 qubits. 
The above p(t) will then be used in calculating various 
entanglement measures. Two qubits in question will be 
either nearest neighbors, that is qubits directly coupled 
by the hamiltonian, or two qubits which are not directly 
coupled, e.g., next nearest neighbors. Typically they will 
be chosen in the middle of the chain with the results 
being largely independent of their precise location. The 
initial pure state will be of two forms. Most of the time 
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we are going to consider random product initial state, 

!#>)> = |X>1®"-®|X>», (12) 
where \x)i is a random state of ith qubit, given by 

\ X )i =cos^e iQ '|0) 4 + 81110,0*11),, (13) 

with Pi, at, 4>i independent random numbers given as 
OLj, = 27r£, Pi = 2tt£ and 4>i = arcsin y/^ where all three £s 
are independent (for each qubit) uniform random num- 
bers on interval [0, 1]. 

IV. ENTANGLEMENT EVOLUTION 
A. Perturbative expansion 

Initially, at time t = 0, our initial state is always of 
product form between two qubits in question and there- 
fore there won't be any entanglement, i.e., C(t = 0) = 0, 
A^ n (i = 0) = and Q(t = 0) = 1. Subsequent evolution 
will entangle two qubits therefore one expects that the 
entanglement will gradually build up. For sufficiently 
short times, one can use perturbation theory to calcu- 
late the reduced density matrix p(t) (|11[) . Taking for 
H nearest-neighbor hamiltonian with an arbitrary two- 
qubit coupling term we get after expanding propa- 
gator exp (—iHt) to the lowest order in time, 

= — 1-5-|* +0(t 2 ), S = (xixi\h (2) \xAXB). (14) 

In the derivation of the above formula we assumed the ini- 
tial product state on the two qubits in question while an 
arbitrary state is allowed on the remaining n — 2 qubits, 
IV'(O)) = \xa) ® \xb) ® \x)n-2- States \xa,b) are sin S le 
qubit states orthogonal to \xa,b) and because only abso- 
lute value of 5 enters their phases do not matter. Similar 
calculation for concurrence © and 9 ^ gives, 

C = 2\S\t + 0(t 2 ) 

9 = l + 4|£|f + 0(t 2 ). (15) 

We can see that the initial speed at which entanglement 
is produced depends only on a single matrix element of 
between the initial product state |xa)<8>|xb) an d the 
corresponding orthogonal product state. Because for all 
quantities the initial time scale depends trivially on the 
value of 1 5 | we will measure time in rescaled dimensionless 
units 

r = t\5\, (16) 

with | S | being the average absolute value of a matrix ele- 
ment, where averaging is done over random initial single 
qubit states \xa,b) (fT3]) . For the isotropic Heisenberg 
model (0 we get \S\ = 1, for the Heisenberg model in a 
staggered field © one gets \S\ « 0.8882 and for the Ising 



model in a tilted magnetic field (J9j) \6\ ~ 0.6168. For two- 
body random matrix model ()10|) we can, instead of aver- 
aging over initial product states, average over ensemble of 
random matrices h^ 2 \ resulting in |<5| = \/tt/^ w 0.4431. 
In all figures showing time dependence of entanglement 
we are going to use dimensionless time r. 

B. Numerical results 
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FIG. 2: Average concurrence C(r) ([2]) in Fig. (a) and aver- 
age 0(r) ([6]) in Fig.(b) for various hamiltonians. Thick full 
curves are for two-body random matrix model (|10[l , thin dot- 
ted curves for the Ising model in a tilted magnetic field ([SJl, 
short dashed curves are for isotropic Heisenberg model (TfJ) 
while long dashed curves are for the Heisenberg model in a 
staggered field ([8| . Averaging is performed over random prod- 
uct initial states (|12[) for n — 16. Times r* when G = 1 are 
0.66 for two-body random matrix model, 0.73 for tilted Ising 
and 0.53 for both Heisenberg models. 

First, we performed numerical simulations of time evo- 
lution, calculating average C(r), A^ n (r) and O(r) for 
different hamiltonians. Averaging has been done over 
random initial product states (|12[) . As one can see in 
Fig. [2] the behavior is overall similar for all studied hamil- 
tonians. For instance, from Fig. [2Jd we can see that for 
times larger than some critical r* the average is be- 
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low 1 which means that the two qubits can not be used 
for distillation any more. Critical t* is for all models 
between 0.5 and 0.75. Similar dependence (not shown) 
is also obtained for the average minimal eigenvalues of 
partially transposed density matrix A^ n , which also be- 
comes positive at roughly the same r* . Concurrence, seen 
in Fig. [2^i, has a similar time dependence. The only dif- 
ference is that for r > r* concurrence is not strictly zero 
but instead decays exponentially with time. This is a con- 
sequence of the fact that even though for t > r* the two 
qubits are on average not entangled any more there are 
still exponentially rare instances (product initial states) 
for which there is still some entanglement present. With 
time the probability of such entangled states decreases 
exponentially In all cases dependence for small times 
agrees with analytical pcrturbative result for concurrence 
in eq.H5land Fig. [2^,, for O in eg . [TBI and Fig.[2j3, and for 
A ml in eq.IHand Fig. H 

Note that r* is a time when the average quantity (like 
or A^f n ) reaches a certain value (1 or 0). It should not 
be confused with the average time f * when (or A^ n 
or C) reaches 1. For each individual initial condition 
time when a state gets separable, i.e. critical t, is of 
course the same as that of C or A^ n . However, becouse 
distributions change with time, time r* is not exactly 
the same as the average time f*. The average times f* 
are for A^f n (as well as for C or negativity) equal to 
f* = 0.72 for two-body random matrix model and tilted 
Ising model and f * = 0.59 for both Hciscnberg models. 
The average times f* for O(r) are on the other hand 
slightly different, f * = 0.64 for two-body random matrix 
model, f* = 0.61 for tilted Ising model and f* = 0.47 
for both Heisenberg models. Compare these f * with r* 
listed in Figs, \2k>, EH and 2] One can observe that f * 
for O and Aj^* n are slightly different, for instance in the 
case of a two-body random matrix model 0.64 vs. 0.72. 
The fact that f * is for O smaller than for A^f n is not a 
contradiction as / < 1/2 does not necessarily mean that 
the entangelement of formation is zero [24| . 

We can see that due to randomness in the initial 
state the evolution of entanglement between two nearest- 
neighbor qubits shows universal-like behavior, that is 
time dependence which is to a large extend independent 
of the precise form of the underlying nearest-neighbor 
hamiltonian which generates evolution. Universality is 
not exact, there are still "signatures" of specific hamilto- 
nian at intermediate times (for instance, compare curves 
for the Heisenberg and Ising model in Fig. overall 
though the dependence is similar. Bell-like shape of en- 
tanglement is a consequence of two competing processes. 
On the one hand, non-separable evolution naturally tends 
to produce entanglement from an initially separable state 
while on the other hand, it will tend to destroy it be- 
cause \i/j(t)) will approach a random state as time grows 
and the reduced density matrix will approach an identity 
matrix having zero entanglement The later process 
of destroying entanglement has two sources: first, ran- 
domness of the initial state is spread out throughout the 



system and second, dynamics itself will tend to produce 
random state. Universality is a consequence of random- 
ness in the initial state, i.e., of its generic separable form. 
For each hamiltonian there are rare specific separable ini- 
tial states for which deviations from the above average 
behavior will be large. 

As we have seen in Fig. [5] two-body random matrix 
model, which is a parameter-free model, describes evolu- 
tion of entanglement sufficiently well also for other sys- 
tems. Therefore, from now on we are going to focus on a 
two-body random matrix model, studying more precisely 
how the entanglement between two qubits evolves with 
time. 



C. Two-body random matrix model 
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FIG. 3: Average Q(t) ([6| for nearest-neighbor qubits (r — 1), 
next nearest (r — 2) and qubits separated by 2 qubits (r = 3). 
For directly coupled qubits, r = 1, states on average can not 
be used for distillation for times larger than r* « 2/3. Chain 
curve (almost overlapping with the full curve for r — 1) is 
® = i+qt'* • ^ U * s f° r two-body random matrix model with 
n = 16 qubits and initial product states (|12|l . In the inset we 
show results for r = 1 and smaller systems. 

In Fig. [3] we again show dependence of ([6]) on time 
for a two-body random matrix model. In addition to 
entanglement between nearest-neighbor qubits (denoted 
by distance r = 1) we also show entanglement for next- 
nearest neighbors (r = 2) and qubits separated by two 
other qubits (r = 3). As one can see, for qubits which 
are not directly coupled by the hamiltonian, i.e., r = 2 
and r = 3 cases, O is always less than 1. This hap- 
pens because the production of entanglement depends 
on higher order terms, e.g., for r — 2 terms of form 
/ii,i+i/ii+i,i+2 are needed, whereas for nearest-neighbors 
a single term /i^j+i is already sufficient to entangle two 
qubits. As a consequence, entanglement production is 
slower the larger is the distance between qubits while 
entanglement destruction due to randomness is approxi- 
mately independent of the distance. In all cases we show 
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data for n — 16 for which finite size effects are already 
small. For instance, the difference between r* for n = 16 
and n = 18 is 0.01 in the case of r = 1. In the figure we 
also plot rational function O(r) = prg!j5 which almost 
perfectly overlaps with the numerics for r = 1. Note 
that for short times this of course agrees with our per- 
turbative result (|15|) . seen as the initial line with slope 4 
in Figs. [2] and [3J One would be tempted to think that 
the dependence beyond this short time, therefore also the 
above rational function, could be explained by higher or- 
der perturbation theory. Unfortunately it is not so. This 
rational dependence can not be explained by higher order 
perturbative calculation. In fact, we do not have any the- 
oretical explanation for this almost perfect fit. Although 
going to perturbations of 2nd order in time will result in 
a rational function with the denominator and numerator 
being polynomials of order 2 in r, the coefficients of poly- 
nomials are wrong. We have numerically checked that 
next perturbative orders also do not improve the situa- 
tion. It therefore seems that the functional dependence 
of (as well as of C and A^ n ) for two-body random 
matrix model is essentially non-perturbative. This is in 
contrast with, for instance, purity or fidelity decay of ini- 
tial pure states in the presence of weak coupling where 
perturbative approaches have been very successful [81.1291]. 
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FIG. 4: Average minimal eigenvalue A m f n of partially trans- 
posed density operator. We show results for nearest-neighbor 
qubits (r = 1), and those for qubits separated by one (r = 2) 
and two (r = 3) qubits. All is for two-body random matrix 
model with n — 16 qubits and initial product states (| 1 2 j) . In 
the inset we show results for r — 1 and smaller systems. 

In Fig. [4] we show numerical results for the average 
^min- Overall, the dependence is very similar as for 0(t), 
the only difference being that the time when be- 

J O mm 

comes positive and the state gets separable is t* r* 0.75. 
In Fig. [5] similar results are shown for concurrence and 
negativity. We can see that negativity and concurrence 
are almost the same for nearest- neighbor qubits. For 
next-nearest neighbor qubits (r = 2) concurrence is this 
time non-zero but small as opposed to which is 



always positive. For qubits further apart, for instance 
r = 3, concurrence is below the level of statistical fluctu- 
ations. 

We have checked that similar results are obtained also 
for other topologies of the coupling between qubits, i.e., 
other than nearest-neighbor. In all situations entangle- 
ment is present for two qubit reduced density matri- 
ces between qubits directly coupled by the hamiltonian, 
whereas entanglement is small or zero for qubits which 
are not directly coupled. In all cases entanglement on av- 
erage disappears after finite time. Difference from sudden 
death of entanglement phenomenon (ll| is that our two 
qubits are coupled and start from an initially separable 
state. In addition, our system is conservative, that is we 
have a finite "environment" . 




FIG. 5: Average concurrence ((2} for nearest-neighbor qubits 
[r = 1) and next nearest-neighbors (r = 2). With dotted line 
we also show negativity, — 2 • N(p) flj), displaying essentially 
the same behavior as concurrence. For two-body random ma- 
trix model with n = 16 qubits and initial product states. 

So far we have always used product initial states, where 
states on individual qubits were independent. As a 
final numerical calculation let us check how the results 
depend on the choice of an initial state. Because results 
are similar for all quantities studied we are going to show 
only 0(t). Besides product initial state (|12p we used 
product initial state with states on all qubits being 
the same. Because two-body random matrix model is 
invariant for single qubit rotations this is equivalent as 
choosing state |0 . . . 0) for the initial state and averaging 
over an ensemble of two-body random matrices. Second 
choice is an initial state which is of product form on the 
two qubits used for entanglement calculation and random 
on the remaining n — 2 qubits, 



|V(0)) = |x)n-2®|XA> 



\XB) 



(17) 



where \xa,b) are random single qubit states and \x)n-2 
is a random state of n — 2 qubits. Numerical results 
are shown in Fig. [6j We can see that, expectedly, en- 
tanglement decays slower for homogeneous initial state, 
|0 . . . 0), while it decays faster for the initial state having 
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a full (nonseparable) random state on n - 2 qubits (fT7f . V. CONCLUSION 

Figuratively speaking one can say that the decay of en- 
tanglement is faster the more randomness there is in the 
initial state. 













... I0...0> 
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i 

FIG. 6: Average O(t) for various choices of initial states. Full 
line is for initial product states (|12|l (same data as in Fig. [3} , 
chain curve is for initial state which is of product form just 
on the two qubits involved (|17[l . while dotted curve is for 
homogeneous initial state |0 . . . 0). All is for two-body ran- 
dom matrix model with n = 16 qubits. Times when average 
6(r) < 1 are 0.42, 0.66 and 0.86. 



We have studied how the entanglement between two 
qubits evolves with time when two qubits are part of a 
larger spin chain. Starting from an initial separable ran- 
dom product state entanglement first increases, reaching 
a maximal value, after which it decays to zero resulting 
in a separable state after finite time. Therefore, starting 
from a generic initial state possessing some randomness, 
two qubit reduced density matrix is on average entangled 
only for finite time and only for qubit pairs which are 
directly coupled by the hamiltonian. Time dependence 
of entanglement is almost independent of the specifics 
of the hamiltonian used for time evolution, being inte- 
grable or chaotic, and is well described by a two-body 
random matrix model. Results can be interpreted also in 
another way: it is hard to generate entanglement regard- 
less of the dynamics if there is some randomness present 
in the initial state. This can be used to explain the lack 
of entanglement in small subsystems for generic initial 
conditions. The author would like to thank T. Prosen 
for reading the manuscript and anonymous referee for 
valuable suggestions and for bringing to our attention 
Refs. [3 and [13 ■ 
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